clear;
data=importdata('adata5.txt');
x=data(:,1);
y1=data(:,2);
[m,~]=size(data);
y2=zeros(m,1);
for i=1:m
    y2(i)=1.0/(1.0+25*x(i)*x(i));
end
themax=max(abs(y2-y1));
disp(themax);
plot(x,y1,'r')
hold on;
plot(x,y2,'b');
